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1. Introduction 

Everyday computational situations are replete with cases where the solution quality does 
not seem to justify the time spent in generating the answer, while a quick even if less 
accurate solution would have been as satisfactory. In the following sections, we describe 
how detection of phase transitions in computational problems may be done, within 
a statistical mechanics framework, thereby avoiding regimes wherein computational 
efforts put in are not commensurate with the solution obtained. This would be of 
utmost importance especially in situations where the problem being solved is an NP 
(nondeterministic polynomial time) hard problem, or where there is a transition in the 
very nature of the problem from an easily computable regime to one that is very hard to 
solve. To be able to actually predict a phase-transition in a computational context (not 
just find it subsequently) and to predict and avoid the "critical slowing down regime" 
has not, to the best of our knowledge, been demonstrated before. 

We discuss this in detail in the context of Anytime Algorithms. This framework 
has been completely explained in this paper itself, so that this work is self-contained 
and complete. We have implemented this method here, yielding very interesting and 
significant results. There is a considerable body of literature dealing with the travelling 
salesman problem and with phase transitions & their use in constraint satisfaction 
problems [I]-[II]. While the use of phase transitions in search problems in general has 
been extensively investigated, there does not seem to have been an attempt in trying 
to use phase transitions in anytime algorithms as a means of monitoring the progress of 
the algorithm or in any other form. We show that a critical exponent can be found that 
relates the time taken to come to a stable solution quality preceding a phase transition, 
to the number of nodes or cities of a travelling salesman problem. We believe this is 
a new and important result that has not been reported in the literature, as far as we 
know. The existence of a critical exponent makes it possible to predict when the onset 
of a computationally hard regime would take place. This is of great importance as it 
has long been unclear how to efficiently run algorithms especially when one has to make 
a trade-off between solution quality and the computational time needed to achieve an 
acceptable, near-optimal solution, given constraints of time. 

The particular problem used - of the Travelling Salesman Problem, is only 
illustrative of the theoretical approach we suggest for detecting phase transitions. 
We show how this relatively simple approach may be used in any other search and 
optimization problem so that computational time may be vastly reduced by halting 
the running of an algorithm when a near-optimal solution is reached prior to a phase 
transition to a computationally hard regime. It would be of interest to see if phase 
transitions could be used in conjunction with anytime algorithms as a means of further 
improving accuracy and speed in searches. 
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2. A brief sketch of anytime algorithms 

Anytime algorithms are algorithms whose execution can be stopped at any time, 
giving a solution to the problem being solved. The quality of the solution improves 
with time [I2l [13]. These algorithms become specially important when we are faced 
with constraints on our resources-which might include, amongst other things, the 
computational resources and the time available to us to solve the problem. Thus, to 
satisfy the constraints put on it, an anytime algorithm trades-off between the solution 
quality and the computational time, for example. Two main types of anytime algorithms 
are contract and interruptible algorithms. Contract algorithms run for a fixed length 
of time, and give a solution only when that interval of time elapses. Interruptible 
algorithms form a larger family - they may be interrupted at any stage to obtain an 
answer though the accuracy and the meaningfulness of the solution would, of course, 
vary. Monitoring an anytime algorithm is of great importance as we can keep track of 
the progress of the algorithm and decide when and at what intervals we could monitor 
or stop the algorithm in order to obtain an optimal answer optimally. 

There are some quantities which play a very basic part in such algorithms. A 
performance profile gives us a measure of the expected quality of the output with 
the execution time. A performance profile is thus typically a probability distribution. 
Those performance profiles which we have tacitly used are actually dynamic performance 
profiles, Pr{Qj\Qi, At), which is the probability of obtaining a solution of quality Qj 
by resuming the algorithm for time interval At when the current solution has quality 
Qi. A Utility function U{Qi,tk) tells us the utility of a solution of quality Qi at time 
tk [IS] . In the work of Hansen and Zilberstein [13] a monitoring policy m is chosen for 
tracking the progress of an algorithm such that it maximises a value function V{fr,tk), 
a cost function Ci being introduced to include the cost of monitoring. As one can only 
estimate the true solution quality at any given time, a feature fr is made the basis for 
estimating the solution quality Qi. Use is made of partial observability functions like 
Pr{Qi\fr,tk) and Pr{fr\Qi,tk) in the value function above to estimate improvements in 
quality. The algorithm is accordingly monitored or allowed to run or halted, as dictated 
by the monitoring policy. Our approach is somewhat different, as explained in the next 
section. 

3. Phase transitions and Anytime Algorithms 

Although phase transitions occur only in the thermodynamic limit, in the limit of in- 
finitely large systems, as was first pointed out by Kramers [H], finite state transitions 
nevertheless exhibit similarities to true phase transitions and hence are very much rel- 
evant to computational situations. Singularities and singular behaviour could be used 
to define universality classes on the basis of common characteristics like, for example, 
swiftly changing correlation lengths between parts of a system at and near the point of 
transition. 
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Phase transitions have appeared in a number of computational contexts and 
paradigms (see, for example |5l El HI], and references therein). The transition from 
a polynomial to exponential search costs, the transition from an underconstrained 
to overconstrained problem, the appearance of transitions in optimization problems, 
automatic planning and models of associative memory amongst others, all indicate the 
widespread prevalence of phase transitions in a computational context. It has also 
been found that on an average several search heuristics have hard problem instances 
concentrated at similar parametric values which points correspond to transitions in 
solubility. 

When one wishes to employ an anytime algorithm for getting a quick, approximate 
solution to a problem, the question that arises is how long is long enough before stopping 
the execution of the algorithm and deciding to accept a solution? Can one decide 
beforehand, on a systematic basis, the length of this run-time? For a problem, the time 
of onset of a transition to a computationally hard regime wherein the solution quality 
does not significantly improve enough to justify running the algorithm for a longer time, 
could act as one criterion for deciding the run-time. Detecting such a transition requires 
one to run the algorithm unmindful of solution quality so that any change in behaviour 
that quantifies the transition can be located. This is what we have done in this paper. 

Another requirement is to describe the particular computational problem being 
considered using some more quantifiable terms. Several authors have addressed 
combinatorial problems through ideas from statistical mechanics; a good review of the 
literature can be found in [11] and references therein. The work of Gent and Walsh [15] 
addresses, in particular, the problem of a phase transition for the Travelling Salesman 
Problem. They identify a transition between soluble and insoluble instances of the 
decision problem (namely, whether or not a tour of some length / or less exists for the 
given TSP), at some critical value of a parameter. 

In physical systems, a transition from fluid to gas states, or from paramagnetic to 
ferromagnetic states for magnetic systems, can be well described using the behaviour 
of the thermodynamical potentials and response functions like susceptibility and order 
parameters like the magnetization (see, for example. Reference [T6]). 
In the following sections, we present a formalism for quantifying phase transitions in 
the computational context. 

The test-problem that we study in this paper is that of obtaining a near-optimal 
solution for a two-dimensional Travelling Salesman Problem (TSP) with A^ nodes or 
cities distributed randomly over a unit square. The premise is that any change in 
behaviour of the computational time of the algorithm would be detected, enabling us 
to quantify any transition to a computationally hard regime. 

The challenge that we face here is how to draw an analogy between thermodynamic 
potentials and response functions for a physical system to quantities in a more abstract 
system such as a computational problem and algorithm. 
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4. Defining tiie formalism in analogy with statistical physics 

In this paper, we consider the two-dimensional Travelhng Salesman Problem as the toy- 
problem to test our formalism on. In keeping with expectations that any phase transition 
in a system should be distinctly observable, we first decide therefore to look for a drastic 
transition in various quantities like the Quality of a solution and the Value function. 
The actual utility of a solution at any time-instant depends upon its quality, as well as 
the rate of improvement shown by the quality and the computational time needed to 
arrive at it. The utility can be expected to increase in direct proportion to both the 
quality of the solution as well as the rate of improvement of quality, while reducing in 
direct proportion to the increasing computational time required. We therefore make 
an ansatz that a Utility function denoted by U{Qi,Qi,tk) can be defined with these 
properties, and which in the simplest instance has a linear dependence on Q and Q: 

U{Qi,Qi,tk) = aiQi + a2Qi- a^t. (1) 

We have taken ai, a2 and 03 to be positive constants. The Quality function Qi can be 
defined in several ways, depending upon the particular problem being considered. For 
the TSP, one obvious choice for Qi would be the reduction in path-length over each 
iteration, since one wishes to complete a tour as optimally as possible. We hence choose 
to define quality by 

g, = (L,-Le)/L„ (2) 

Li being the initial path length and Lc the current path length. The value function 
(mentioned earlier in Section 2) used is taken as a sum of the expected values of the 
utility function, and is defined by 

V{U.tu) = Y.^Pr{Qj\U,tu,M)U{Q,,Q,,tk + M). (3) 

In the toy model being discussed, we will take the feature fk to be the quality Qk-, 
so that V{Qk,tk), defines in quantitative terms the value of a solution of quality Qk 
at time tk in obtaining other solutions of quality Qj in a time step following tk- The 
utility U{Qj, Qj, tk + At) when weighted by the conditional probability of obtaining that 
solution of quality Qj in time step At after starting at time tk with quality Qk, and 
summed over all Qj, clearly gives a measure of the value of the solution of quality Qk- 
Hence the definition (3) of the value function. 

The procedure adopted is as follows. 2-Opt heuristics is employed to solve each 
problem instance. To calculate and generate a performance profile, the algorithm is 
run for about 500 times. To then solve a particular problem, the performance profile 
generated earlier is used on completion of each iteration of the algorithm to calculate 
quantities like the Value function. The initial tour is generated by means of a nearest- 
neighbour algorithm. The steps of the initial tour construction are: 

(i) Select a random city, 
(ii) Find the nearest unvisited city and go there. 
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(iii) If any unvisited cities are left, repeat the previous step, 
(iv) Return to the first city. 

Tour improvement is done through a 2-Opt heuristic; i.e., following initial tour construc- 
tion, in subsequent iterations, the algorithm selects a small segment of the sequence and 
reverses it. We use a simple 2-bond move which reverses the sequence of cities for a 
chosen segment to obtain a trial tour. In each iteration of this 2-opt heuristic, we accept 
a tour if it is shorter in length than the previous one, rejecting it otherwise. Solution 
quality therefore keeps improving or remains unchanged, as long as time permits. On 
completion of each iteration of our algorithm, we use the performance profile generated 
earlier to calculate the value function and the cost inflection function described below. 
Data from 500 samples are averaged for obtaining each point in the data sets used to 
draw our inferences. 
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Figure 1. Frequency of occurence of different Quality values for an algorithm for a Trav- 
elling Salesman Problem for number of nodes A^ = 30, 50 and 70, over 10® instances. The 
distribution obeyed is eqn.(4), with a, /3 & 7 values being approximately given by: 72.57, 
43.67, 255.45, respectively for iV = 30; k 242.72, 46.31, 425.31, respectively for N = 50; and 
317.69, 51.55, 617.56, respectively for N = 70. 



The frequency with which a certain quality Q is achieved is shown in Figure 1 for 
a Travelling Salesman tour with the number of nodes N = 30, 50, and 70. The quality 
is found to obey a distribution of the form 

P(Q) = aexp(/3Q-7Q3), (4) 

a, (3 and 7 being iV-dependent, positive parameters. 

The question thus arises - how do we locate or predict a phase transition in a problem 
such that it can be of practical import and used? To deal with this, we define a function 
K, which we call the Cost Inflection function, the suddenly changing behaviour of which 
would indicate the occurence of a phase transition. Since we want as good a solution 
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as soon as possible, with ininimal time elapse, we could think of K as the value of 
the solution corrected to account for the time expended in arriving at it. We make 
the ansatz that the Cost Inflection function K{Qi,ti) is related to the Value function 
V{Qi,t) and to solution quality through the relation: 

K{Qi,t) = V{Qi,t)-tQ„ (5) 

where the value V{Qi,t) is defined through equation (3), Qi being the quality function 
of the currentlv available solution. 
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Figure 2. Plot of the Value function V 
as a function of time t, for different cases 
of Travelling Salesman Problems with 
number of nodes N = 100, 250, 300 and 350. 



Figure 3. Three dimensional plot of Cost- 
Inflection, Quality and Value Functions, 
for N = 150, 250 and 350. 



Figures 2 and 3 show the dependence between Cost Inflection, Value and Quality 
functions. As can be clearly seen, the Cost Inflection function shows a marked change 
in behaviour after a particular point. We identify K as an analogue to the Helmholtz 
free energy A. 

Recall that the Value function is a weighted sum of the Utility function. The 
Value function's temporal evolution (Figure 2) reminds us of the time evolution of 
the most probable value Xm of the number of particles (corresponding to the chemi- 
cal composition) for a system in the stochastic theory of adiabatic explosion [17]. In 
that thermodynamic system, X^ evolves with time till some critical time tc after which 
new solution-branches appear corresponding to other probability states. In a system 
undergoing combustion, this is the temporary situation where the molecules can be dif- 
ferentiated into a part for which combustion has not yet taken place and a part for 
which combustion has ended. In our problem, the evolution of V{fi,t) is an evolution 
of the system to higher quality solutions, differentiated by a point of inflection at a time 
t = tc where the system slows down heralding the onset of critical slowing down. 
The appearance of a symmetric spike in the V profile plotted as a function of t and 
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K, and a sharp inflection of tlie curve into anotlier plane, as sliown in Figure 4, are 
symptomatic of a major transition in quality, and is accompanied by a transition to a 
different computational regime. 



N=250(K> 




Figure 4. The second derivative of the 
Value Function with respect to time, V 
plotted as a function of K and time t, the 
Quality of the solution, Q, being depicted 
by the colour; number of nodes N = 250. 
Note the sharp inflection corresponding 
to a phase transition. 
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Figure 5. The Quality Function plotted as a 
function of time for three values of N {N = 
100, 200, 250), and an overlaying plot of an 
Arrhenius function, to show the broad 
qualitative similarity. 



In Figure 5 we have plotted the improvement of solution quality with time for the 
TSP defined above, for different values of N. We find that the behaviour of Q can be 
approximated by means of an Arrhenius equation 

Qapprox = beyip{-c/t), (6) 

b and c being some real constants. The suffix approx has been added to Q to stress 
that this is but a tool in our gedankenexperiment to help understand a rather abstract 
system in terms more familiar to us. This form of Qapprox is almost exactly like the 
Arrhenius law for the rate constant w{T) in a chemical reaction in a combustive process 
w{T) = y4exp(^), E being the activation energy, R the gas constant, T the temper- 
ature, and A the pre-exponential or frequency factor that can be approximated to be 
constant. Just as in the thermodynamic situation the rate constant increases rather 
steeply with temperature, eventually plateauing out to a constant value, in our case 
the solution quality too increases quickly with time before approaching the maximum 
achievable value for the problem. 
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Figure 6. Representative K — t — K — Q plots obtained for travelling salesman problems with 
AT = 30, AT = 100, N = 2m k N = 300. K is the Cost Inflection Function, t the time, K 
corresponds to the second derivative of K with respect to time. Q, the Quality function, is 
represented by the colour. 



As in thermodynamics, in our computational context too one can define appropriate 
response functions. As time elapses, the system responds by a perceptible change in its 
quality. This is reflected as a decrease in the cost inflection function K after a critical 
value (see Figure 6). This motivates us to define a quantity Cat, for which we coin the 
term efficacy of the algorithm: 

CN = \tk\, (7) 

where derivatives are with respect to time. The cost inflection function K is calculated 
from equation (5). A plot of the efficacy function. Cat, as a function of time t for 
different A^ values (Figure 7 a-b) shows an interesting observation - Cj^ tends to show 
almost random behaviour upto a point where a transition to a stable regime takes place; 
thereafter, the efficacy function shows a strict t dependence 

C-N stable OC t. (8) 
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Figure 7. Efficacy Cn plotted as a function of time t (a) for different number of 
nodes N, all showing a transition to a stable regime where C^ stable oc t^; (b) for A^ = 175 
at left, and N = 300 at right: representative plots showing the regime of Cat = CNstabie 
to correspond to the high quality metastable & stable states. 



This stable area, as can be seen in Figure 7b, corresponds to regimes where the 
quality Q has reached its maximal value, or is reaching it, and encompasses, therefore, 
both stable and metastable states. The efficacy as defined above therefore makes a very 
good response function in identifying a transition to a near-optimal set of regimes. 
In the thermodynamic context, the specific heat C^ for a fluid is, as we know, defined 
by 



G^, 



-T 



gj2 



(9) 
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where A denotes the Helmholtz free energy. The specific heat behaviour as a function of 
temperature is distinct, depending upon whether the temperature is below or above a 
critical temperature Tc for the thermodynamic system. For example, for the well-known 
superfluid or A transition in "^He the specific heat dependence can be approximated by 
C ~ A log(T — Tx) + B, for temperatures above the critical A-transition temperature 
T > Tx, and C ^ A' log(Tx ~ T) + B' for temperatures below the critical temperature 
T < Tx, with the prefactors A, B, A' and B' being temperature - independent 
constants [18] . 

In our case, the response function Cjy is such that its behaviour changes abruptly at 
the point corresponding to a critical value of elapsed time for the system. In other 
words, looking at Cjy, we should clearly be able to distinguish between separate phases 
or regimes. 

Approach to a phase transition or bifurcation point in physical systems is typically 
characterized by a critical slowing down. We therefore expect that our system should 
likewise show similar behaviour. In the Travelling Salesman Problem under consid- 
eration, critical slowing down would correspond to a regime wherein solution quality 
remains unchanged over a comparatively long period of time, while other measures of 
algorithm efficiency, like the cost inflection function K, or the value function V, tend 
to show a change in behaviour. An investigation of our data indicates that the onset of 
critical slowing down in the algorithm can be quantified and detected. Locating tc, the 
time of onset of this critical slowine down, is quite straightforward. 




N 

Figure 8. Plot indicating dependence of t^ time of onset of critical slowing down, on 
number of nodes A'^. The best fit for the log- log plot is the equation tc = a^N^, where 
aO = 0.00047 ± 0.00027, z = 2.073 ± 0.099. 

It should be noted that there may be more than one regime where the solution 
quality of the algorithm lingers at one value for an extended period of time. Such states 
may essentially be thought of as metastable states for the system, as has been mentioned 
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in our discussion of the efficacy function above, tc, however, is distinguished from such 
metastable states by the clearly discernable change in behaviour of functions like K, V 
or V. We find that a plot of tc, the time corresponding to the onset of critical slowing 
down, versus A^, the number of cities or nodes in a 2-D Travelling Salesman Problem, 
(Figure 8), scales like: 

tc ~ iV^ 2 ^ 2.07. (10) 

This critical exponent is reminiscent of the dynamical scaling exponent z, in physical 
systems [lEl HH |19] . We note that our system is intrinsically non-autonomous. Figure 
6 shows representative plots of K plotted against K and time, with the Quality of the 
solution being represented by the colour. This is yet another way of looking at how 
there is a change in the evolution of the Cost Inflection function with time. 
The presence of large fluctuations can be noted in the plots shown in Figure 6 and, for 
example, in Figure 4. Apart from the fact that such fluctuations would be expected as 
one approaches a point of bifurcation. Figure 7a suggests too that for a given value of 
A^, the system can have more than one, close "energy state", so that the system tends 
to fall into either one or more of these metastable and stable states during the course 
of its evolution. Generic theoretical treatments of fluctuations near bifurcations that 
correspond to the appearance of new stable states have appeared in the literature (for 
example. Reference [20]). However, an analysis of the fluctuations observed in our sys- 
tem and the apparent absence of correlation of fluctuation size to system size, is outside 
the scope of this present work, and will not be dealt with here. 

Differentiating eqn.(5) twice with respect to time, 

tQ + 2Q = V-K. (11) 

Since K — V is very small and near zero for most times after the critical point 
corresponding to the time t^ beyond which we can justifiably expect that no further 
temporal evolution of the rate at which the Value and Cost Inflection functions change 
with time takes place, we can write 

tQ + 2Q^ 0. (12) 

Integrating this between the limits t = ta and tb, tf, > ta, we get 

Qb-Qa = H^~h^ (13) 

ta tb 

k being some constant of integration. If we now identify ta as the time corresponding 
to critical slowing down, tc, and take tb to be much larger, so that tb — )■ cxd, we are left 
with 

where Qc is the value of the quality function at the onset of critical slowing down 
and Qmax is the value when the algorithm is run for a much longer time. It will be 
appreciated that equations (13) and (14) are valid only for time t > tc- 
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As we have seen (Figure 5, equation (6)), the variation of Q with time t can be 
very broadly approximated by an Arrhenius-hke functional dependence on time Q = 
6exp(— c/t). At large times, this can be approximated to Q ^ b{l — c/t), so that 

We are therefore able to give a value to the constant k of equation (14), 

k ^ be, (16) 

so that the solution quality at t^ is readily predictable: 

^;c ^max C IZq 

^6(l-c/g, (17) 

Smce ycf^rnax ^^ "■ 

An order parameter rUg for this time-evolving system could be a measure of quality 
defined by 

mg = Q,-k{l/t-l/t,). (18) 

For t > tc, nig would necessarily be greater than equal to Qc, i.e., nig — Qc>0. For any 
time t < tc, niq would be less than Qc, so that niq — Qc < 0. 

Another issue of interest in travelling salesman problems is obtaining an estimate 
of the optimal tour length. The optimal tour distance L^ for a d-dimensional TSP is 
a function of d and the number of cities N. In the large N limit, this optimal tour 
distance is given by 



ao being a constant. Percus and Martin have investigated the finite size scaling of the 
optimal tour length of this Euclidean TSP with randomly distributed cities [22]. For 
the two-dimensional case, they found that Lo can be written in the form 

iVV.(l + V(8Ar) + .„) =°°"-°-°^^^/^ + -'- P") 

with ao = 0.7120 [22]. In an earlier paper reporting the simulated annealing simulations 
of a travelling salesman problem [23] , Lee and Choi obtained the optimal tour length to 
be given by 

Lo{N) = aVN + (3, (21) 

a = a{N —7- oo), a{N) = Lo/vN, getting an upper bound value of a = 0.7211, 
(3 = 0.604. 

It would also be interesting to know how close to optimality are the solutions that 
we have obtained from our algorithm. Keeping in mind that what we seek to do in 
this paper is to describe a technique for deciding when to halt a computation while 
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approaching a solution that is approaching optimahty, rather than obtaining optimum 
solutions, we find from our data that for our best solutions, 

Lo{N) ^ 0.7408ViV + 1.1732 (22) 

holds, i.e., a = 0.7408, and /3 = 1.1732. Figure 9 shows this in a plot of the length of 
the oDtimal solution L„ versus l/\/N. 
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Figure 9. Plot of optimal length Lq vs. 1/\/N . The plot obeys equation (22), Lo{N) 
0.7408\/iV + 1.1732. 



The general trend of our solutions, therefore, are consistent with that in the 
literature. Improved values for a can more easily be obtained when the system size 
is much larger, and when computations are run longer. 

5. Discussion and Conclusions 

From the above treatment, one can see that it would be possible to look out for drastic 
transitions in search costs and solution quality while monitoring algorithms. In the 
case of Anytime Algorithms, for example, we could now effectively identify transitions 
from contract to interruptible algorithms. Identifying the transition and the associated 
behaviour, required us of course, to run our algorithm for a long period of time, 
irrespective of solution quality. It should be noted that the quality measure could be 
defined in several ways, and we have chosen but one possible definition. Our choice of 
obtaining an initial configuration using a greedy assignment of nearest neighbour jumps 
was made because it works well for a random 2D Travelling Salesman Problem, which 
was the toy problem we chose for testing our ideas. This need not, of course, be used 
for other problems. 

Our method does not set out to explicitly find a global minimum corresponding 
to some particular stable state, but rather, to locate the existence of both stable 
and metastable states. The intention is to be able to find and predict regimes where 
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computational time spent is not commensurate with improvement in solution quality 
achieved. 

We have defined a response function, somewhat akin to a specific heat of a fiuid, 
whose behaviour handily identifies metastable and stable states corresponding to high 
solution qualities. This efficacy Cn is a clear indicator of a transition from lesser quality 
unstable solutions to the higher quality stable solutions, showing a t-dependence for the 
latter, for the TSP problem we have considered. 

We have identified a dynamical critical exponent z ^ 2.07 relating the onset of 
critical slowing down place to the number of nodes in a travelling salesman problem. 
This is a very invaluable result: we can now know when a transition in quality will occur 
and have a handle on how long an algorithm need be run. It will be recalled that all 
possible rearrangements in a 2-opt algorithm are of order A^^. The coincidence of this 
with our critical exponent of 2.07 underlines the fact that our approach is robust and is 
consistent with what would be expected from a more traditional approach. It is crucial, 
however, to note that it is not the enumeration of A^^ 2-opt states that we have called 
a dynamical critical phenomenon, but rather the onset of critical slowing down, since 
time is indeed a control parameter with all our functions (quality, cost, utility, response 
function, etc.) varying dynamically in time. We expect that similar universal exponents 
would be valid for other examples in this class of problems. 

It becomes meaningful therefore to conduct similar investigations on other problems 
and look for similar critical exponents. Further work on this is underway and will be 
reported elsewhere. 

While we have implemented our formalism and demonstrated its working using a 2-opt 
algorithm for a travelling salesman problem, this generic form may be suitably modified 
for use on any other problem using any other algorithm. This is the real advantage of 
our work - it enables one to set the time parameters under which an anytime algorithm 
may be used as a contract anytime algorithm, i.e., with the time of its running being 
fixed in advance; we are also able to quantify an abstract quantity like solution quality 
in a generic fashion by means of a response function which clearly indicates the onset 
of a phase transition in solution quality and hardness. 
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